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Abstract.  This  paper  describes  CHEM2D-H20,  a  new  pa¬ 
rameterization  of  JLO  photochemical  production  and  loss 
based  on  the  CHEM2D  photochemical-transport  model  of 
the  middle  atmosphere.  This  parameterization  accounts  for 
the  altitude,  latitude,  and  seasonal  variations  in  the  photo¬ 
chemical  sources  and  sinks  of  water  vapor  over  the  pressure 
region  from  100-0.001  hPa  (~16-90km  altitude).  A  series 
of  free-running  NOGAPS-ALPHA  forecast  model  simula¬ 
tions  offers  a  preliminary  assessment  of  CHEM2D-H20  per¬ 
formance  over  the  June  2007  period.  Results  indicate  that 
the  CHEM2D-H20  parameterization  improves  global  10- 
day  forecasts  of  upper  mesospheric  water  vapor  compared  to 
forecasts  using  an  existing  one-dimensional  (altitude  only) 
parameterization.  Most  of  the  improvement  is  seen  at  high 
winter  latitudes  where  the  one -dimensional  parameterization 
specifies  photolytic  PLO  loss  year  round  despite  the  lack 
of  sunlight  in  winter.  The  new  CHEM2D-H20  parameter¬ 
ization  should  provide  a  better  representation  of  the  down- 
welling  of  dry  mesospheric  air  into  the  stratospheric  polar 
vortex  in  operational  analyses  that  do  not  assimilate  middle 
atmospheric  FLO  measurements. 


1  Introduction 

Although  the  middle  atmosphere  (15- 100  km  altitude)  is  ex¬ 
tremely  dry  when  compared  to  the  troposphere,  detailed 
knowledge  of  the  water  vapor  distribution  in  this  region  is 
important  for  a  number  of  reasons.  For  example,  water  vapor 
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abundance  controls  the  availability  of  odd  hydrogen  species 
for  catalytic  ozone  loss.  In  addition,  the  emission  of  long¬ 
wave  (terrestrial)  radiation  to  space  by  water  vapor  is  an  im¬ 
portant  cooling  process  in  the  middle  atmosphere.  The  rel¬ 
atively  long  photochemical  lifetime  of  middle  atmospheric 
water  vapor  also  makes  it  a  useful  tracer  for  studying  the  dy¬ 
namics  of  this  region.  Finally,  the  abundance  of  middle  at¬ 
mospheric  water  vapor  is  an  important  factor  controlling  the 
formation  of  both  polar  stratospheric  clouds  in  winter  and 
polar  mesospheric  (or  noctilucent)  clouds  near  the  summer 
mesopause,  the  former  being  important  for  heterogeneous 
ozone  loss.  In  more  practical  terms,  water  vapor  is  a  fun¬ 
damental  prognostic  variable  in  the  dynamical  cores  of  most 
numerical  weather  prediction  (NWP)  models.  For  these  rea¬ 
sons,  NWP  and  data  assimilation  (DA)  systems  whose  top 
levels  extend  into  the  upper  stratosphere  and  mesosphere  re¬ 
quire  an  accurate  description  of  the  photochemical  sources 
and  sinks  of  water  vapor. 

In  general,  operational  requirements  for  timely  forecasts 
prevent  NWP  systems  from  performing  fully  coupled  pho¬ 
tochemical  model  calculations  because  they  are  too  compu¬ 
tationally  intensive.  One  solution  to  this  problem  is  employ 
parameterizations  based  on  linearization  of  the  relevant  pho¬ 
tochemical  processes  (see,  e.g.  Cariolle  and  Deque,  1986). 
For  a  review  of  the  development  of  linearized  photochem¬ 
istry  parameterizations,  see  McCormack  et  al.  (2006). 

This  paper  describes  a  new  linearized  gas-phase  water  va¬ 
por  photochemistry  parameterization  that  is  based  on  the 
CHEM2D  zonally  averaged  photochemical-transport  model 
of  the  middle  atmosphere.  CHEM2D  has  been  successfully 
used  to  develop  fast,  accurate  parameterizations  of  strato¬ 
spheric  ozone  photochemistry  (McCormack  et  al.,  2004, 
2006).  MacKenzie  and  Harwood  (2004)  implemented  this 
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J.  P.  McCormack  et  al.:  Parameterized  PPO  photochemistry 


(a)  CHEM2D  CH4  loss  rates 


(b)  CHEM2D  H20  loss  rates 


Fig.  1.  Vertical  profiles  of  the  dominant  CHEM2D  model  loss  rates  (molecules  cm  3  s  at  5°  N  on  15  June  for  (a)  CH4  and  (b)  FPO. 


type  of  water  vapor  photochemistry  parameterization  in  a 
general  circulation  model  to  study  trends  in  middle  atmo¬ 
spheric  humidity.  Here  we  apply  a  similar  approach  to  pro¬ 
vide  a  suitable  representation  of  middle  atmospheric  water 
vapor  for  operational  NWP/DA  systems. 

This  new  water  vapor  photochemistry  parameterization, 
designated  CHEM2D-H20,  has  been  recently  implemented 
in  the  high-altitude  version  of  the  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS-ALPHA).  Here 
we  present  a  description  of  CHEM2D-H20  as  well  as  re¬ 
sults  from  a  series  of  NOGAPS-ALPHA  forecast  model 
simulations  designed  to  evaluate  CHEM2D-H20  and  com¬ 
pare  its  performance  to  a  simpler  one-dimensional  (1-D)  wa¬ 
ter  vapor  photochemistry  scheme  currently  used  in  opera¬ 
tional  NWP/DA  systems.  The  ultimate  goal  of  this  work 
is  to  provide  global  simulations  of  mesospheric  water  va¬ 
por  accurate  enough  to  identify  the  physical  processes  gov¬ 
erning  polar  mesospheric  cloud  (PMC)  formation  near  the 
summer  mesopause.  Section  2  gives  a  general  overview 
of  the  CHEM2D  model  and  middle  atmospheric  water  va¬ 
por  photochemistry.  Section  3  describes  the  implementa¬ 
tion  of  CHEM2D-H20  in  NOGAPS-ALPHA.  Section  4  ex¬ 
amines  middle  atmospheric  water  vapor  fields  from  a  se¬ 
ries  of  NOGAPS-ALPHA  forecast  model  simulations  to  as¬ 
sess  the  performance  of  CHEM2D-H20.  Section  5  summa¬ 
rizes  these  results  and  outlines  future  research  applications 
of  CHEM2D-H20. 


2  H2O  photochemistry  in  the  CHEM2D  model 

CHEM2D  is  a  zonally  averaged  (2-D)  global  model  that 
features  a  fully  self-consistent  treatment  of  radiative,  pho¬ 
tochemical,  and  dynamical  processes  of  the  middle  atmo¬ 
sphere  (see,  e.g.,  McCormack  et  al.,  2006;  McCormack  et 
al.,  2007).  The  model  photochemistry  accounts  for  reactions 
among  54  different  species  using  reaction  rates  from  Sander 


et  al.  (2003).  CHEM2D  extends  from  pole  to  pole  with  grid 
points  spaced  every  4.8°  in  latitude;  the  vertical  domain  con¬ 
sists  of  88  pressure  levels  from  the  surface  to  p= 6  x  10-5  hPa 
(~  1 1 6  km)  spaced  every  ~  1 .3  km. 

CHEM2D  computes  the  solar  UV  heating  above  40  km  ex¬ 
plicitly  using  the  same  spectral  distribution  of  solar  UV  and 
absorption  cross  sections  as  in  the  model  photolysis  calcula¬ 
tions.  Details  of  the  photolysis  calculations  can  be  found  in 
Summers  et  al.  (1997)  and  references  therein.  Modeled 

O2  absorption  in  the  Schumann-Runge  bands  is  based  on 
the  parameterization  of  Minschwaner  et  al.  (1993).  H2O  ab¬ 
sorption  is  treated  as  in  Siskind  et  al.  (1994).  O2  and  H2O 
absorption  at  Lyman-a  is  based  on  the  method  of  Lewis 
et  al.  (1983).  The  effects  of  CH4  absorption  at  Lyman-a 
on  H2O  are  not  considered.  Solar  cycle  variations  in  irra- 
diance  are  specified  from  1200-8000  Abased  on  the  mea¬ 
surements  of  Lean  et  al.  (1997).  The  CHEM2D  photoly¬ 
sis  rates  reported  here  are  for  solar  minimum  conditions. 
The  Lyman-a  irradiance  at  solar  minimum  was  taken  to  be 
3xl0n  photons  cm-1  s  . 

Water  vapor  (H2O)  in  the  middle  atmosphere  is  produced 
directly  through  oxidation  of  stratospheric  methane  (CH4) 
by  the  hydroxyl  radical  (OH+CH4  — >  H2O+CH3)  and  in¬ 
directly  through  a  series  of  reactions  involving  the  methyl 
radical  (CH3).  The  net  effect  is  that  approximately  two 
H2O  molecules  are  produced  for  each  CH4  molecule  lost  in 
the  stratosphere  (e.g.,  LeTexier  et  al.,  1988,  and  references 
therein).  Consequently,  the  total  number  density  of  hydrogen 
throughout  much  of  the  stratosphere  Q=2[CH4 1+|  H2O]  is  a 
constant  (neglecting  the  relatively  small  amount  of  molecu¬ 
lar  hydrogen,  H2).  This  relationship  can  be  used  to  express 
stratospheric  H2O  production  in  terms  of  H2O  abundance, 
which  is  advantageous  since  it  eliminates  the  need  for  a  prog¬ 
nostic  CH4  variable. 

Figure  la  plots  vertical  profiles  of  individual  CHEM2D 
CH4  loss  rates  (molecules  cm-3  s_1)  at  5°  N  on  15  June.  Re¬ 
actions  (5)  and  (6)  listed  in  Fig.  la  are  the  key  reactions 
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(a)  CHEM2D  rCH4  Jun 


(b)  CHEM2D  rCH4  Dec 


96 

80 


Latitude 


(c)  CHEM2D  tH20  Jun  (d)  CHEM2D  rH20  Dec 


Fig.  2.  Pressure-latitude  plots  of  CHEM2D  CH4  photochemical  lifetimes  (in  days)  for  (a)  15  June  and  (b)  15  December,  and  of  CHEM2D 
H2O  lifetimes  for  (c)  15  June  and  (d)  15  December.  Shading  in  (a)  and  (b)  indicates  region  where  CH4  loss  is  primarily  through  photolysis. 
Contours  are  drawn  at  3,  5,  10,  30,  50,  100,  300,  500,  and  1  000  days. 


responsible  for  stratospheric  EDO  production.  Above  the 
0. 1  hPa  level,  CH4  loss  occurs  through  photolysis  and  does 
not  produce  EDO  directly. 

As  a  result  total  hydrogen  Q  is  not  conserved  in  the 
mesosphere.  Middle  atmospheric  H2O  is  destroyed  pri¬ 
marily  through  photolysis  in  the  mesosphere  by  solar  UV 
and  Lyman-a  radiation,  yielding  a  number  of  different  by¬ 
products  depending  on  the  energy  of  the  incoming  radia¬ 
tion  (Steif  et  al.,  1975).  EDO  is  lost  in  the  stratosphere 
through  reaction  with  atomic  oxygen.  Figure  lb  plots  in¬ 
dividual  CHEM2D  vertical  profiles  of  EDO  loss  rates  at  5°  N 
on  15  June.  The  loss  rates  for  the  photolysis  Reactions  1- 
3  in  Fig.  lb  are  based  on  estimated  quantum  yields  of  0.78, 
0.10,  and  0.12,  respectively.  These  values  of  the  quantum 
yields  were  determined  from  a  synthesis  of  available  labora¬ 
tory  measurements  (Steif  et  al.,  1975;  Harich  et  al.,  2000). 

Photochemical  loss  rates  are  commonly  used  to  infer  an 
effective  photochemical  lifetime,  which  is  a  convenient  way 
to  quantify  the  relevant  time  scales  for  photochemistry  rela¬ 
tive  to  other  physical  processes  such  as  advection.  Figure  2 
plots  CHEM2D  photochemical  lifetimes  of  CH4  and  EbO 
(tch4  and  rnzo  respectively)  for  June  and  December  condi¬ 
tions  throughout  the  middle  atmosphere.  The  lifetimes  are 
computed  from  the  sum  of  the  individual  loss  rates  as 


[CH4]  [H20] 

TCH4  =  ^6  T  ’  riFO  =  4  , 

2-,i=\  2^i=\  L-i 


(1) 


where  [CH4]  and  [H2O]  denote  CHEM2D  CH4  and  EDO 
abundances  (molecules  cm-3)  and  L;  is  the  loss  rate 
(molecules  cm~3  s“ 1 )  for  the  individual  reactions  in  Fig.  1 . 

Figures  2a  and  b  show  that  the  shortest  CH4  lifetimes 
(<2  days)  are  in  the  summer  mesopause  region  where  pho¬ 
tolysis  is  the  dominant  loss  mechanism.  A  secondary  mini¬ 
mum  in  tcH4  is  seen  near  the  summer  polar  stratopause  due 
to  rapid  CH4  loss  via  reaction  with  chlorine  (Reaction  4  in 
Fig.  1). 

The  shaded  regions  in  Fig.  2  indicate  where  CH4  loss  is 
dominated  by  photolysis,  a  process  that  does  not  directly  pro¬ 
duce  H2O.  The  unshaded  regions  of  Fig.  2a  and  b  indicate 
where  CH4  loss  produces  EDO  in  the  stratosphere  via  Reac¬ 
tions  5  and  6  as  in  Fig  la. 

Values  of  rcm  for  levels  below  0. 1  hPa  exceed  50  days,  in¬ 
dicating  that  stratospheric  EDO  production  is  slow  compared 
to  typical  transport  timescales.  As  a  result,  including  the 
effects  of  this  process  in  typical  NWP/DA  systems  issuing 
5-  to  10-day  forecasts  is  not  crucial  provided  that  accurate 
stratospheric  humidity  data  are  being  regularly  assimilated 
globally.  In  the  absence  of  such  data,  or  when  conducting 
longer  free-running  model  simulations  for  seasonal  predic¬ 
tion  or  climate  simulations,  the  contribution  of  CH4  oxida¬ 
tion  to  EDO  production  becomes  more  significant. 

Fig.  2c  and  d  show  that  values  of  rmo  in  the  summer 
hemisphere  are  less  than  5  days  above  the  0.01  hPa  level, 
indicating  that  EbO  photolysis  is  an  important  effect  for 
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(a)  CHEM2D  (P-L)H20  Jun  (b)  CHEM2D  (P-I_)H20  Dec 


Fig.  3.  Pressure-latitude  plots  of  net  CHEM2D  middle  atmosphere  water  vapor  photochemical  tendency  (P— L)h20  f°r  (a)  15  June  and 
(b)  15  December  in  parts  per  million  (volume)  per  month.  Contours  are  drawn  at  ±0.1,  ±0.3,  ±1,  ±3,  and  ±10ppmv  per  month.  Dashed 
contours  denote  negative  values. 


r)r 


accurate  medium-range  forecasts  in  this  region.  Values  of 
t"H20  near  the  mesopause  increase  rapidly  poleward  of  50° 
latitude  in  winter  as  the  amount  of  sunlight  diminishes. 

To  highlight  the  latitude,  altitude,  and  seasonal  depen¬ 
dences  of  H2O  photochemistry  in  the  middle  atmosphere. 
Fig.  3  plots  the  net  CHEM2D  photochemical  tendency  for 
water  vapor,  (P— L)h2Cu  in  parts  per  million  by  volume 
(ppmv)  per  month  as  a  function  of  latitude  and  pressure  for 
15  June  and  15  December,  where  P  is  the  net  production  rate 
and  L  is  the  net  loss  rate.  In  the  stratosphere,  (P—L) H20 
is  weakly  positive  as  a  result  of  relatively  slow  FbO  pro¬ 
duction  through  CH4  oxidation  in  the  absence  of  water  va¬ 
por  photolysis  at  these  altitudes  (see  Figs.  1  and  2).  In  the 
mesosphere,  (P— L)h20  is  mostly  negative  due  to  photolytic 
loss,  with  the  largest  loss  rates  occurring  near  the  summer 
mesopause  region.  The  positive  values  of  (P— T)h20  in 
the  winter  hemisphere  near  80  km  poleward  of  50°  latitude 
are  the  result  of  enhanced  FbO  production  via  the  reaction 
OH+HO2  — >  H2O+O2.  This  enhancement  is  due  to  pole- 
ward  transport  of  odd  hydrogen  from  sunlit  latitudes  into  po¬ 
lar  night,  where  it  becomes  long-lived  in  the  upper  meso¬ 
sphere  (Brasseur  and  Solomon,  1986). 

The  values  of  tch4-  Dqo,  and  (P— L)h20  presented  in 
this  section,  which  vary  with  latitude,  altitude,  and  month, 
serve  as  the  basis  for  the  new  CHEM2D-FbO  parameteri¬ 
zation.  CHEM2D-H2O  differs  from  the  water  vapor  photo¬ 
chemistry  parameterization  currently  used  in  the  European 
Centre  for  Medium  Range  Weather  Forecasts  (ECMWF)  In¬ 
tegrated  Forecast  system  (IFS)  (Untch  and  Simmons,  1999; 
ECMWF,  2006;  Feist  et  al.,  2007),  which  parameterizes  pro¬ 
duction  via  stratospheric  CH4  oxidation  and  loss  via  meso¬ 
spheric  photolysis  as  a  function  of  altitude  only,  neglecting 
possible  latitude  and  seasonal  dependences.  This  method  ex¬ 
presses  the  water  vapor  photochemical  tendency  as 


—  =k\(rQ-r) -k2r  (2) 

at 

where  r  is  the  local  water  vapor  mixing  ratio  and  tq  is  the 
equivalent  total  hydrogen  mixing  ratio. 

The  coefficients  k\  and  ki  are  determined  from  analytical 
fits  to  quoted  values  of  tch4  in  the  stratosphere  and  TH20 
in  the  mesosphere,  respectively,  at  various  altitudes  from 
Brasseur  and  Solomon  (1986). 

The  coefficients  k\  and  kj  vary  with  altitude  and  are  con¬ 
stant  with  latitude  and  season.  It  is  assumed  that  i'q  has  a 
constant  value  of  6.8  ppmv  based  on  the  results  of  Randel 
et  al.  (1998).  Figure  4  plots  CHEM2D  values  of  rg  show¬ 
ing  this  assumption  holds  for  most  of  the  stratosphere  except 
winter  polar  regions  where  downward  transport  of  dry  meso¬ 
spheric  air  occurs  (LeTexier  et  al.,  1988;  Randel  et  al.,  1998). 
Thus  the  ECMWF  scheme  provides  a  reasonable  1-D  pho¬ 
tochemical  constraint  on  global  stratospheric  FbO,  but  may 
overestimate  stratospheric  FbO  production  in  regions  where 
conservation  of  i-q  breaks  down. 

To  compare  this  ID  approach  to  the  CHEM2D-H20  pa¬ 
rameterization,  Fig.  5  plots  the  vertical  profile  of  the  com¬ 
bined  photochemical  lifetime,  (Aq+fo)-1  from  the  ECMWF 
scheme  (ECMWF,  2006)  along  with  CHEM2D  values  of 
rcH4  and  th20  over  the  equator  for  each  month  of  the  year. 
Only  CHEM2D  values  of  rcH4  below  the  0. 1  hPa  level  are 
plotted  in  Fig.  5,  since  it  is  in  this  region  where  CH4  loss 
leads  to  FRO  production.  Also  plotted  in  Fig.  5  are  values  of 
the  combined  CHEM2D  lifetime 


rCH4  +  TH20 

for  each  month,  which  are  determined  from  the  plotted  values 
of  rCH4  and  rH20  in  Fig.  5. 

The  quantity  r*  provides  a  single,  concise  description  of 
the  time  scale  for  FEO  photochemistry  throughout  the  middle 
atmosphere  that  varies  with  latitude  and  season  analogous  to 
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Fig.  4.  Pressure-latitude  plots  of  total  hydrogen  rg  (in  ppmv)  from  the  CHEM2D  model  for  (a)  15  June  and  (b)  15  December.  Contours  are 
drawn  every  1  ppmv. 
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Fig.  5.  Comparison  of  combined  photochemical  lifetime  (days)  in 
the  ECMWF  water  vapor  photochemistry  parameterization  (black 
curve)  with  CHEM2D  lifetimes  rjj2o  (purple  curves),  rcu4  (gold 
curves),  r*  (red  curves)  for  all  12  months  over  the  equator.  The 
CHEM2D  rcH4  profile  is  only  plotted  at  levels  where  CHq  oxida¬ 
tion  is  the  dominant  loss  process  (c.f.,  Fig.  1  and  unshaded  region 
in  Fig.  2). 

the  1-D  profile  of  (k\  +/.'2)~ 1  in  the  ECMWF  scheme.  The 
overlapping  curves  in  Fig.  5  indicate  that  there  is  virtually  no 
seasonal  dependence  in  the  FFO  production  and  loss  terms 
at  low  latitudes,  in  agreement  with  Figs.  2  and  3.  As  a  result, 
CHEM2D  values  of  r*  for  all  12  months  closely  match  the 
ID  scheme’s  k j-1  and  kk 1  profile. 

At  higher  latitudes,  however,  seasonal  variations  in  th20 
and  tcH4  are  both  quite  large.  As  a  result,  r*  also  exhibits 
large  seasonal  variations  at  higher  latitudes.  To  illustrate 
this  point.  Figure  6  compares  values  of  (ki+k2)~l  from  the 
ECMWF  scheme  with  values  of  r*  at  70°  N  for  March,  June, 
September,  and  December. 


Fig.  6.  Comparison  of  combined  photochemical  lifetime  (days)  in 
the  ECMWF  water  vapor  photochemistry  parameterization  (black 
curve)  with  values  of  the  CHEM2D  effective  lifetime  r*  (red 
curves)  at  70°  N  on  the  15th  day  of  March,  June,  September,  and 
December. 

The  CHEM2D  model  results  presented  in  this  section 
demonstrate  that  the  net  middle  atmospheric  EDO  photo¬ 
chemical  tendency  exhibits  pronounced  latitude  and  seasonal 
variations  that  should  be  accounted  for  in  NWP/DA  systems 
extending  into  the  mesosphere.  The  CHEM2D-H20  param¬ 
eterization  described  in  the  following  section  accounts  for 
these  variations. 

3  CHEM2D-H20  in  NOGAPS-ALPHA 

This  section  describes  the  CHEM2D-H20  parameteriza¬ 
tion  and  its  implementation  in  NOGAPS-ALPHA.  First,  we 
briefly  describe  the  NOGAPS-ALPHA  forecast  model  and 
data  assimilation  components.  We  then  provide  a  detailed 
description  of  the  CHEM2D-H20  scheme  as  it  is  currently 
used  in  NOGAPS-ALPHA. 
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3.1  NOGAPS-ALPHA  description 

The  high-altitude  NOGAPS-ALPHA  NWP/DA  system 
(Hoppel  et  al.,  2008)  is  capable  of  assimilating  middle  atmo¬ 
sphere  constituent  measurements  such  as  H2O  profiles  ob¬ 
tained  from  the  NASA  EOS  Aura  Microwave  Limb  Sounder 
(MLS)  (Eckermann  et  al.,  2008,  referred  to  hereafter  as 
E2008).  The  68-level  (L68)  NOGAPS-ALPHA  global  spec¬ 
tral  forecast  model  (GSFM)  initializes  and  advects  spe¬ 
cific  humidity  q  at  all  levels  from  the  surface  to  its  top  at 
5xlO_4hPa.  The  q  fields  are  initialized  using  traditional 
meteorological  analyses  from  the  surface  up  to  the  200  hPa 
level.  Above  this  level,  the  humidity  fields  can  be  initial¬ 
ized  using  either  assimilated  water  vapor  measurements  for  a 
specific  date  or  a  zonal  monthly  mean  climatology  based  on 
measurements  from  the  Upper  Atmospheric  Research  Satel¬ 
lite  (UARS)  Halogen  Occultation  (HALOE)  and  MLS  instru¬ 
ments  (Randel  et  al.,  1998;  Gmol.i  and  Russell,  2005),  de¬ 
pending  on  the  application. 

Parameterizations  for  moist  physics  in  NOGAPS-ALPHA 
are  identical  to  those  used  in  the  operational  version  of  NO¬ 
GAPS  (Hogan  and  Rosmond,  1991).  These  include  shallow 
cumulus  mixing  (Tiedtke,  1984),  deep  cumulus  convection 
(Peng  et  al.,  2004),  and  convective,  stratiform,  and  bound¬ 
ary  layer  cloud  formation  and  precipitation  (Slingo,  1987; 
Teixeira  and  Hogan,  2002).  For  a  more  complete  descrip¬ 
tion  of  NOGAPS-ALPHA  model  physics,  see  McCormack 
et  al.  (2006)  and  E2008.  Since  they  are  designed  primarily 
for  tropospheric  applications,  the  model’s  moist  physics  rou¬ 
tines  are  only  employed  from  the  surface  up  to  the  50hPa 
level  (~20  km).  In  the  stratosphere  and  mesosphere,  pa¬ 
rameterized  water  vapor  photochemical  production  and  loss 
constrain  the  NOGAPS-ALPHA  q  fields.  Without  this  pho¬ 
tochemical  constraint,  upwelling  of  stratospheric  air  into 
the  mesosphere  would  produce  unrealistically  high  values 
of  mesospheric  humidity  over  forecast  periods  of  2-3  days. 
This  would  pose  a  problem  for  accurate  forecasts  of  meso¬ 
spheric  humidity,  especially  in  regions  where  humidity  mea¬ 
surements  are  not  assimilated. 

All  forecast  model  results  in  this  study  employ  triangular 
truncation  of  the  first  79  wavenumbers  (T79),  with  a  hori¬ 
zontal  grid  spacing  of  1.5°  in  longitude  and  latitude  on  the 
quadratic  Gaussian  grid.  The  model  uses  the  non-orographic 
gravity  wave  drag  parameterization  of  Garcia  et  al.  (2007) 
with  the  same  settings  as  in  E2008.  The  T79L68  model  is  ini¬ 
tialized  with  analyzed  winds,  temperature,  and  constituents 
(e.g.,  O3  and  H2O)  produced  by  the  NOGAPS-ALPHA  DA 
component,  which  is  described  in  detail  by  Hoppel  et  al. 
(2008)  and  E2008.  From  the  surface  to  the  mid-stratosphere, 
these  NOGAPS-ALPHA  analyses  are  based  on  assimilation 
of  conventional  meteorological  data  sets  used  by  the  opera¬ 
tional  T239L30  system. 

For  the  June  2007  period  studied  here,  NOGAPS-ALPHA 
assimilates  Aura  MLS  temperature,  O3,  and  H2O  profiles 
(Froidevaux  et  al.,  2006),  as  well  as  temperatures  from  the 


Sounding  of  the  Atmosphere  Using  Broadband  Radiome- 
try  (SABER)  instrument  on  the  Thermosphere  Ionosphere 
Mesosphere  Energetics  and  Dynamics  (TIMED)  satellite 
(Kutepov  et  al.,  2006),  up  to  the  0.002  hPa  level,  (E2008). 
Above  this  level,  model  fields  are  initialized  with  climatol¬ 
ogy  as  in  McCormack  et  al.  (2006).  The  altitude  resolution  of 
the  MLS  water  vapor  retrievals  decreases  significantly  with 
increasing  altitude  (Lambert  et  al.,  2007).  A  Gaussian  aver¬ 
aging  kernel  in  the  vertical  direction  with  a  full-width  at  half 
maximum  value  of  2  km  is  applied  to  the  assimilated  MLS 
water  vapor  profiles.  This  was  done  to  constrain  the  analysis 
to  the  low  resolution  MLS  profiles  rather  than  risk  introduc¬ 
ing  any  spurious  vertical  structure  in  the  H2O  fields  that  may 
arise  in  data-poor  regions.  The  background  error  variance  for 
water  vapor  in  the  upper  atmosphere  was  set  large  enough  to 
ensure  that  the  analysis  was  highly  constrained  to  the  obser¬ 
vations,  producing  a  zonal  mean  distribution  that  is  nearly 
identical  to  the  MLS  data.  For  further  details  of  the  assimila¬ 
tion  system,  see  Hoppel  et  al.  (2008). 

A  comprehensive  evaluation  of  the  NOGAPS  middle  at¬ 
mospheric  data  assimilation  system,  including  comparisons 
with  independent  observations,  is  currently  underway.  The 
purpose  of  this  paper  is  to  introduce  the  CHEM2D-H20  pho¬ 
tochemistry  parameterization  and  evaluate  its  performance 
based  on  NOGAPS-ALPHA  forecast  model  simulations  of 
middle  atmospheric  water  vapor.  The  following  section  de¬ 
scribes  the  implementation  of  CHEM2D-H20  in  NOGAPS- 
ALPHA. 

3.2  The  CHEM2D-H20  parameterization 

The  CHEM2D-H20  parameterization  in  NOGAPS-ALPHA 
expresses  the  local  time  rate  of  change  of  the  H2O  mixing 
ratio  r  as  the  difference  between  the  zonally  averaged  pro¬ 
duction  and  loss  rates  computed  with  the  CHEM2D  model 
(see  Fig.  3): 

dr 

7T"  =  (p  ~  L) H20-  (4) 

at 

Although  specific  humidity  q  is  the  NOGAPS-ALPHA 
model’s  prognostic  variable,  for  the  sake  of  consistency  with 
the  photochemical  parameterization  in  Eq.  (2)  we  will  use 
volume  mixing  ratio  r  as  the  moisture  variable  in  the  follow¬ 
ing  discussion.  Model  mixing  ratio  r  is  related  to  q  through 
the  relation 


where  M,j  and  Mw  are  the  molecular  weights  of  dry  air  and 
H2O,  respectively. 

We  assume  that  (P— L)h20  is  primarily  a  function  of  r. 
In  the  mesosphere,  this  assumption  is  justified  because  the 
loss  rate  due  to  photolysis  is  directly  proportional  to  the  local 
H2O  mixing  ratio.  In  the  stratosphere,  where  the  production 
rates  via  methane  oxidation  depend  on  the  local  CH4  mixing 
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Fig.  7.  Pressure-latitude  plots  of  the  effective  photochemical  lifetime  r*  (days)  used  in  the  CHEM2D-H20  parameterization  for  (a)  15  June 
and  (b)  15  December. 


ratio,  this  assumption  holds  where  the  local  FPO  mixing  ra¬ 
tio  can  be  expressed  in  terms  of  the  CH4  mixing  ratio  through 
the  approximate  conservation  of  i-q  (see  Fig.  4).  Where  the 
assumption  of  rg  conservation  breaks  down,  i.e.,  in  the  win¬ 
ter  polar  stratosphere,  CHEM2D  values  of  rgm  exceed  100 
days.  As  shown  below,  accounting  for  these  latitude  and 
seasonal  variations  in  FbO  production  and  loss  ensures  that 
the  CHEM2D-H20  parameterization  only  introduces  photo¬ 
chemical  tendencies  in  regions  where  its  underlying  assump¬ 
tions  are  valid. 

The  function  (P—L) H20M  can  be  approximated  using 
a  first-order  Taylor  series  expansion  about  a  reference  state 
such  that 


d r(X,  <p,  p,  t) 
3t 


(P~L)°u 


H20 


9(-P  -  L) H2P 

dr 


(r  -  r°)  (6) 


where  X  is  longitude,  (p  is  latitude,  p  is  pressure,  and  “o” 
denotes  the  reference  state. 

We  proceed  by  assuming  that  Eq.  (6)  yields  an  equilibrium 
(reference)  state  r°  that  is  the  net  balance  between  photo¬ 
chemical  production  and  loss,  such  that  any  deviations  from 
this  state  can  be  treated  as  small  perturbations  about  that 
reference  state.  This  follows  the  method  used  in  linearized 
ozone  photochemistry  schemes,  as  reviewed  by  McCormack 
et  al.  (2006)  and  Lahoz  et  al.  (2007). 

Applying  the  expansion  Eq.  (6)  to  the  ECMWF  relation 
Eq.  (2),  it  is  straightforward  to  show  that 


{P-L)°H2o  =k1(rQ-r°)-k2r0,  (7) 

d(P-L)°n 

- —*°\0  =  -(kl+k2)  (8) 


illustrating  that  this  approach  is  mathematically  consistent 
with  the  1-D  parameterization  in  Eq.  (2). 

As  discussed  in  Sect.  2  and  illustrated  in  Figs.  5  and  6,  the 
CHEM2D-H20  analog  of  Eq.  (8)  is 


d(P-L) 


h2o 


dr 


0  =  -(**)" 


(9) 


where  r*  represents  the  combined  time  scales  for  FhO  pro¬ 
duction  via  CH4  oxidation  in  the  stratosphere  and  FhO  pho¬ 
tolysis  in  the  mesosphere.  CHEM2D  model  values  of  r*  are 
tabulated  as  a  function  of  latitude,  pressure  and  season  and 
then  interpolated  in  space  and  time  to  the  NOGAPS -ALPHA 
forecast  model  grid.  Figure  7  plots  the  latitude  and  altitude 
dependence  of  r*  for  June  and  December  conditions. 

The  analytical  ECMWF  expression  Eq.  (7)  assumes  con¬ 
stant  total  hydrogen  rg=6.8ppmv  at  all  altitudes.  As  illus¬ 
trated  in  Fig.  4,  this  assumption  breaks  down  in  the  upper 
stratosphere  and  mesosphere.  In  addition,  CH4  oxidation  be¬ 
comes  a  negligible  source  of  H2O  in  the  mesosphere.  Thus, 
for  CHEM2D-H20  we  instead  derive  (P— L)^i0  values  as 
a  function  of  latitude,  pressure  and  season  directly  from  the 
CHEM2D  model,  examples  of  which  were  shown  previously 
in  Fig.  3.  This  approach  enables  us  to  compute  equilibrium 
rates  without  any  explicit  reference  to  total  hydrogen,  and 
thus  to  extend  these  rates  into  the  mesosphere  where  total  hy¬ 
drogen  conservation  breaks  down  and  is  dominated  by  pho- 
tolytic  loss. 

Adopting  an  approach  similar  to  that  for  the  CHEM2D 
ozone  photochemistry  parameterization,  or  CHEM2D-OPP 
(McCormack  et  al.,  2006),  the  specific  humidity  photochem¬ 
ical  tendency  in  NOGAPS-ALPHA  is  applied  by  first  defin¬ 
ing  a  photochemical  steady  state  value  for  the  water  vapor 
mixing  ratio 

rss  =  r°  +  (P  -  L)° r*  (10) 


so  that  the  H2O  mixing  ratio  tendency  can  be  expressed  as 
dr  _  —  (r  -  rss ) 
dt  r* 


(11) 


The  updated  mixing  ratio  value  is  computed  using  a  standard 
backward-Euler  solution  of  the  form 


r(t  +  At)  =  r(t)  +  [rss  -  r(t)] 


Al 


l  +  #f 


(12) 


and  then  converted  to  specific  humidity  using  Eq.  (5). 


www.  atmos-chem-phys. net/8/75 1 9/2008/ 


Atmos.  Chem.  Phys.,  8,  7519-7532,  2008 


7526 


J.  P.  McCormack  et  al.:  Parameterized  PPO  photochemistry 


Fig.  8.  Pressure-latitude  plots  of  the  zonal  mean  reference  state  H2O  mixing  ratio  r°  (in  ppmv)  for  June  based  on  (a)  the  combined  UARS 
HALOE/MLS  climatology  used  in  EXP1  and  (b)  the  Aura  MLS  climatology  used  in  EXP2. 


Table  1.  FLO  photochemistry  parameterizations  used  in  the 
NOGAPS-ALPHA  forecast  model  simulations. 


Experiment  Name 

Description 

EXP1 

CHEM2D-H20  with  r° 

based  on  UARS  HALOE/MLS 

EXP2 

CHEM2D-H20  with  r° 

based  on  Aura  MLS 

EXP3 

ECMWF 

In  theory,  the  reference  state  mixing  ratio  r°  should  cor¬ 
respond  to  the  CHEM2D  model  mixing  ratio  at  photochem¬ 
ical  equilibrium.  In  practice,  values  of  r°(ip,  p,  t)  are  often 
specified  using  an  observation-based  climatology  of  middle 
atmospheric  water  vapor.  Based  on  earlier  experience  with 
linearized  ozone  photochemistry  parameterizations  (McCor¬ 
mack  et  al.,  2006;  Geer  et  al.,  2007;  Coy  et  al.,  2007),  use 
of  an  observation-based  reference  state  ensures  that  the  lin¬ 
earized  photochemical  tendency  terms  will  not  produce  large 
biases  between  the  modeled  and  assimilated  constituent  val¬ 
ues,  which  can  negatively  impact  forecast  skill.  In  this  study, 
we  perform  a  similar  evaluation  of  the  CHEM2D-H20  pa¬ 
rameterization’s  sensitivity  to  the  assumed  background  state 
using  two  different  r°  distributions.  The  first  is  based  on 
monthly  zonal  mean  MLS/HALOE  climatology  (Randel  et 
al.,  1998;  Grool.i  and  Russell,  2005)  between  100-0.01  hPa 
combined  with  CHEM2D  model  values  above  the  0.01  hPa 
level.  The  second  combines  June  zonal  mean  Aura  MLS 
H2O  mixing  ratios  averaged  over  the  years  2005-2008  up  to 
the  0.002  hPa  level  with  CHEM2D  model  values  above  this 
level.  Solar  UV  fluxes  in  the  CHEM2D  model  were  set  to 
solar  minimum  levels  following  Lean  et  al.  (1997)  to  match 
June  2007  conditions. 


Figure  8  compares  these  two  r°  distributions.  The  Aura 
MLS  H2O  mixing  ratios  in  Fig.  8b  are  generally  higher  than 
the  UARS-based  values  in  Fig.  8a  throughout  the  Northern 
Hemisphere  extratropical  upper  mesosphere.  The  results  pre¬ 
sented  in  the  following  section  demonstrate  how  these  dif¬ 
ferences  can  affect  10-day  forecasts  of  middle  atmospheric 

h2o. 

4  Results 

Table  1  lists  three  sets  of  NOGAPS-ALPHA  forecast  model 
simulations,  designated  EXP1,  EXP2,  and  EXP3,  used  to 
evaluate  the  new  CHEM2D-H20  parameterization.  EXP1 
employs  CHEM2D-H20  with  the  UARS-based  r°  for  the 
month  of  June  (see  Fig.  8a).  EXP2  employs  CHEM2D- 
H20  with  the  Aura  MLS-based  r°  for  June  (Fig.  8b).  EXP3 
uses  the  ID  water  vapor  photochemistry  parameterization  (2) 
currently  used  in  the  operational  ECMWF  IFS.  These  three 
sets  of  simulations  each  consist  of  five  10-day  forecasts  ini¬ 
tialized  on  5,  10,  15,  20,  and  25  June  2007,  encompassing 
the  period  when  PMC’s  were  observed  in  the  Arctic  region 
(E2008).  The  EXP1,  EXP2,  and  EXP3  simulations  all  use 
the  same  sets  of  initial  conditions. 

We  first  examine  the  performance  of  CHEM2D-H20  us¬ 
ing  the  UARS-based  reference  state  (EXP1).  Figure  9a  plots 
the  water  vapor  mixing  ratio  initial  conditions  for  00:00  UTC 
on  5  June  2007,  which  are  based  on  the  assimilation  of 
Aura  MLS  Version  2.2  (V2.2)  H2O  profiles  between  100- 
0.002  hPa  (~  16-90  km)  combined  with  CHEM2D  model 
values  of  the  H2O  above  the  0.002  hPa  level  (see  E2008  for 
details).  Note  that  EXP1,  EXP2,  and  EXP3  simulations  all 
use  the  same  sets  of  initial  conditions. 
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Fig.  9.  Pressure-latitude  plots  of  (a)  NOGAPS-ALPHA  analyzed  zonal  mean  FPO  mixing  ratios  for  00:00  UTC  5  June  2007;  (b)  zonal  mean 
forecast  model  H2O  mixing  ratios  at  day  10  of  the  EXP1  model  simulation  initialized  00:00  UTC  5  June  and  valid  00:00  UTC  15  June  2007 
(c)  zonal  mean  NOGAPS-ALPHA  analyzed  H2O  at  00:00  UTC  on  15  June  2007;  (d)  differences  in  H2O  mixing  ratios  between  the  10-day 
forecast  (F)  in  (b)  and  the  analyzed  (A)  values  in  (c).  Solid  and  dashed  contours  denote  positive  and  negative  differences,  respectively. 


Figure  9b  plots  the  model  H2O  mixing  ratios  on  day  10  of 
the  EXP1  simulation  valid  15  June  at  00:00  UTC.  A  com¬ 
parison  of  Fig.  9b  with  Fig.  9a  shows  that  the  major  differ¬ 
ences  in  NOGAPS-ALPHA  zonal  mean  H2O  after  10  days 
are:  (1)  a  decrease  at  all  latitudes  above  the  0.01  hPa  level 
(~82km);  (2)  an  increase  near  0.1  hPa  (~65km)  poleward 
of  50°  N;  (3)  a  decrease  over  high  southern  latitudes  between 
1 .0-0. 1  hPa.  The  high  latitude  changes  are  consistent  with 
upwelling  (downwelling)  of  moist  (dry)  air  in  the  summer 
(winter)  polar  mesosphere.  The  broad  decrease  above  the 
0.01  hPa  level  is  the  result  of  the  parameterized  photolysis  of 

h2o. 

Figure  9c  plots  the  zonal  mean  NOGAPS-ALPHA  an¬ 
alyzed  H2O  mixing  ratios  valid  on  00:00  UTC  15  June 
2007.  This  and  all  following  plots  of  analyzed  H2O  ex¬ 
tend  to  0.0025  hPa,  which  is  the  closest  vertical  level  to 
the  0.002  hPa  upper  limit  for  scientifically  useful  MLS  H2O 
measurements  (Froidevaux  et  al.,  2006).  Since  NOGAPS- 
ALPHA  analyzed  H2O  is  completely  constrained  by  the 
MLS  measurements,  the  analyzed  zonal  mean  values  are  es¬ 
sentially  the  same  as  a  zonal  mean  estimated  directly  from 
the  daily  MLS  data.  Comparison  of  the  analyzed  H2O  with 
the  10-day  fields  from  EXP1  (Fig.  9b)  shows  large  differ¬ 
ences  poleward  of  30°  N  above  the  0.01  hPa  level,  where  the 
analyzed  values  are  much  higher  than  the  modeled  values. 
Figure  9d  plots  the  differences  between  the  zonal  mean  ana¬ 


lyzed  H2O  and  10-day  forecast  H2O.  After  10  days  the  EXP1 
model  simulation  underestimates  the  zonal  mean  H2O  mix¬ 
ing  ratio  throughout  the  northern  extratropical  upper  meso¬ 
sphere  by  more  than  2  ppmv.  This  is  the  same  region  where 
the  UARS-based  and  Aura  MLS-based  r°  distributions  differ 
substantially  (Fig.  8),  and  where  photolytic  loss  dominates 
(Fig.  3). 

As  mentioned  in  the  previous  section,  the  performance  of 
the  CHEM2D-H20  parameterization  can  be  affected  by  the 
choice  of  the  background  reference  state  r°  because  it  re¬ 
laxes  the  forecast  model  H2O  mixing  ratios  in  the  upper  lev¬ 
els  toward  the  reference  state  value  r° .  In  practice,  the  water 
vapor  distribution  used  to  represent  r°  should  be  chosen  so 
as  to  avoid  any  systematic  bias  between  the  prognostic  hu¬ 
midity  variable  and  the  assimilated  humidity  fields  used  to 
initialize  and  update  the  forecast  system  (see,  e.g.  Coy  et  al„ 
2007).  Systematically  low  values  of  r°,  such  as  those  seen 
in  Fig.  8a,  can  lead  to  an  overestimate  of  the  H2O  loss  in  the 
northern  summer  upper  mesosphere.  This  may  help  explain 
why  the  10-day  EXP1  forecast  underestimates  the  H2O  mix¬ 
ing  ratio  compared  to  the  analyses  (Fig.  9d)  in  the  northern 
extratropical  upper  mesosphere. 

To  investigate  this  possibility,  we  isolate  the  effect  of  pa¬ 
rameterized  water  vapor  photochemistry  by  taking  the  differ¬ 
ence  between  the  model  mixing  ratio  r  and  a  passive  humid¬ 
ity  tracer  rpass,  which  uses  the  same  initial  conditions  as  r 
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(a)  EXP  1  Ar  (ppmv)  day  =  1 0. 
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(b)  EXP  2  Ar  (ppmv)  day  =  10. 


Fig.  10.  Zonal  mean  values  of  Ar,  the  difference  between 
NOGAPS- ALPHA  model  prognostic  H2O  and  passive  (i.e.,  no  pho¬ 
tochemistry)  H2O  on  day  10  of  the  (a)  EXP1,  (b)  EXP2,  and  (c) 
EXP3  simulations  initialized  00:00  UTC  5  June  2007. 

but  is  not  subject  to  parameterized  photochemistry  and  moist 
physics.  In  the  troposphere  and  lowermost  stratosphere,  the 
lack  of  moist  physics  in  the  passive  humidity  field  produces 
huge  differences  between  r  and  rpass  (not  shown).  In  the  up¬ 
per  stratosphere  and  mesosphere,  differences  between  r  and 
rpass  are  due  primarily  to  the  effect  of  the  CHEM2D-H20 
parameterized  photochemistry  acting  on  r. 


Figure  10  plots  zonal  mean  values  of  the  difference 
A  r—r  —  rpass  throughout  the  middle  atmosphere  on  day  10  of 
EXP1,  EXP2,  and  EXP3  simulations  initialized  00:00  UTC 
5  June  2007.  In  all  three  simulations,  the  effects  of  the  pa¬ 
rameterized  photochemistry  on  the  model  H2O  after  10  days 
are  relegated  to  the  upper  mesosphere.  Specifically,  values 
of  Ar  are  negative  throughout  the  Northern  Hemisphere  be¬ 
tween  0.03-0.001  hPa,  consistent  with  H2O  loss  via  photoly¬ 
sis  (e.g..  Fig.  5)  that  peaks  at  high  summer  latitudes.  The  ex¬ 
tension  of  negative  values  of  Ar  across  the  equator  and  into 
the  winter  hemisphere  is  a  consequence  of  photochemic ally- 
processed  air  with  lower  mixing  ratio  r  being  transported 
from  the  northern  summer  mesosphere  to  the  southern  win¬ 
ter  hemisphere  by  the  mean  meridional  circulation.  Since  the 
timescale  for  H2O  production  via  CH4  oxidation  is  >50  days 
in  most  of  the  stratosphere  (see  Fig.  5),  values  of  Ar  after 
10  days  in  Fig.  10  are  quite  small  in  magnitude  throughout 
the  stratosphere. 

Comparing  the  amount  of  photochemical  loss  among  the 
three  different  simulations  in  Fig.  10,  we  find  that  the  great¬ 
est  amount  of  loss  is  produced  in  EXP1  (Fig.  10a).  The 
amount  of  loss  produced  by  CHEM2D-H20  is  reduced  in 
EXP2  when  the  UARS-based  climatological  values  of  r° 
are  replaced  with  the  Aura  MLS  monthly  zonal  mean  H2O 
(Fig.  10b),  particularly  between  the  0.01  and  0.001  pres¬ 
sure  levels.  Both  EXP2  and  EXP3  results  show  similar 
amounts  of  loss  in  the  upper  mesosphere  at  high  northern  lati¬ 
tudes.  This  indicates  that  both  CHEM2D-H20  and  ECMWF 
parameterizations  produce  comparable  results  in  the  region 
where  peak  H2O  photolysis  occurs.  Figure  10c  also  indicates 
that  the  ECMWF  parameterization  used  in  EXP3  produces 
loss  at  high  southern  latitudes  where  no  sunlight  is  present, 
unlike  the  CHEM2D-H20  results  in  EXP1  and  EXP2. 

To  further  characterize  the  performance  of  CHEM2D- 
H20  in  NOGAPS-ALPHA,  we  compare  10-day  forecast 
model  H2O  fields  with  analyzed  H2O  using  five  sets  of  sim¬ 
ulations  initialized  00:00  UTC  on  5,  10,  15,  20,  and  25  June 
2007.  Figure  11a  plots  vertical  profiles  of  the  global  area- 
weighted  mean  difference  between  analyzed  and  forecast 
H2O  mixing  ratios  (denoted  “A-F”)  from  all  five  sets.  We 
find  that  the  CHEM2D-H20  parameterization  with  UARS- 
based  r°  values  in  EXP1  produces  the  largest  mean  values 
of  A-F  above  the  0. 1  hPa  level,  consistent  with  the  excessive 
H2O  loss  seen  in  Fig.  10.  In  contrast,  CHEM2D-H20  with 
Aura  MLS-based  r°  values  produces  the  smallest  mean  A-F 
values.  This  result  holds  true  when  we  examine  mean  A- 
F  values  over  separate  latitude  bands  between  90°  S-90°  N 
(Fig.  llb-f).  Figure  11  also  shows  that  the  ECMWF  scheme 
in  EXP3  generally  produces  larger  values  of  A-F  than  the 
CHEM2D-H20  scheme  in  EXP2,  particularly  in  the  south¬ 
ern  extratropics  (Fig.  lib  and  c).  This  is  due  to  the  fact 
that  the  ECMWF  scheme  specifies  photolytic  loss  of  H2O  at 
the  high  southern  (winter)  latitudes  despite  the  lack  of  sun¬ 
light.  At  high  northern  (summer  latitudes)  CHEM2D-H20 
and  ECMWF  H2O  results  are  comparable. 
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Fig.  11.  Mean  differences  between  NOGAPS-ALPHA  analyzed  H2O  and  10-day  forecast  H2O  (A-F)  based  on  forecast  model  simulations 
initialized  00:00  UTC  on  June  5,  10,  15,  20,  and  25,  2007  using  CHEM2D-H20  with  UARS-based  reference  state  (EXP1),  CHEM2D- 
H20  with  Aura  MLS-based  reference  state  (EXP2),  and  the  ECMWF  parameterization  (EXP3).  Results  are  presented  for  (a)  the  global 
area-weighted  average  A-F,  (b)  60°  S-90°  S.  (c)  30°  S-60°  S,  (d)  30°  S-30°  N,  (e)  30°  N-60°  N,  and  (f)  60°  N-90°  N. 


The  preliminary  evaluation  of  forecast  model  results  for 
the  June  2007  period  indicates  that  all  of  the  10-day  forecasts 
tend  to  underestimate  the  FbO  mixing  ratios  in  the  northern 
extratropics  (see  Fig.  lie)  compared  to  the  analyzed  MLS 
values  above  the  0.01  hPa  level,  regardless  of  the  details  of 
the  photochemistry  parameterization.  It  should  be  noted  that 
the  Aura  MLS  accuracy  limits  (±0.5  ppmv)  and  broad  (12- 
16  km)  vertical  resolution  in  the  upper  mesopheric  FLO  re¬ 
trievals  (Lambert  et  al.,  2007)  may  complicate  the  evalua¬ 
tion  of  forecast  skill  in  this  region.  It  should  also  be  noted 
that  the  forecast  model  transport  in  the  upper  mesosphere  is 
highly  sensitive  to  the  details  of  the  parameterized  gravity 
wave  drag,  and  that  these  details  remain  poorly  constrained 
by  observations.  Because  the  effective  FLO  lifetime  is  com¬ 
parable  to  transport  timescales  in  the  0.1-0.01  hPa  region, 
deficiencies  in  model  transport  may  contribute  to  some  of  the 
disagreement  between  the  observations  and  the  10-day  fore¬ 
cast  model  simulations.  The  effects  of  model  transport  in 
general,  and  gravity  wave  drag  in  particular,  on  NOGAPS- 
ALPHA  middle  atmospheric  H2O  simulations  are  currently 
under  investigation. 


Finally,  to  demonstrate  the  effect  of  the  CHEM2D-H20 
and  ECMWF  photochemistry  parameterizations  on  strato¬ 
spheric  water  vapor,  where  the  relevant  photochemical  time 
scales  are  much  longer  than  in  the  mesosphere,  the  EXP1 
and  EXP3  simulations  were  extended  out  to  day  90.  Fig¬ 
ure  12  compares  zonal  mean  A r  values  from  EXP1  and 
EXP3  between  10-0. 1  hPa  on  day  90.  In  general,  we  find  that 
both  schemes  produce  small  increases  in  water  vapor  (0.25- 
0.5  ppmv)  in  the  tropical  upper  stratosphere.  A  notable  result 
in  Fig.  12  is  that  the  CHEM2D-H20  scheme  in  EXP1  leads 
to  relatively  more  photochemical  loss  of  H2O  in  the  meso¬ 
sphere  compared  to  the  ECMWF  scheme  in  EXP3  and  thus 
a  drier  upper  stratosphere  over  the  South  Pole  due  to  down¬ 
ward  transport  of  mesospheric  air  at  high  winter  latitudes. 
This  result  suggests  that  using  a  H2O  photochemistry  pa¬ 
rameterization  with  latitude  and  seasonal  dependences  may 
reduce  the  moist  bias  in  mesospheric  H2O  seen  in  ECMWF 
analyses  (Feist  et  al.,  2007),  and  thereby  improve  its  repre¬ 
sentation  of  dry  mesospheric  air  descending  into  the  polar 
vortex  during  winter  and  spring.  To  address  this  issue  further, 
we  plan  to  compare  NOGAPS-ALPHA  water  vapor  analyses 
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(a)  EXP  1  Ar  (ppmv)  day  =  90  (b)  EXP  3  Ar  (ppmv)  day  =  90 
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Fig.  12.  Zonal  mean  values  of  the  difference  between  NOGAPS-ALPHA  prognostic  H2O  and  passive  FPO  mixing  ratio,  Ar,  on  day  90  of 
(a)  EXP1  and  (b)  EXP3  free-running  NOGAPS-ALPHA  forecast  simulations  initialized  00:00  UTC  5  June  2007. 


from  the  full  assimilations  system  using  both  the  CHEM2D- 
H20  and  ECMWF  parameterizations  over  the  course  of  the 
entire  season. 

5  Conclusions 

We  have  described  and  tested  the  new  CHEM2D-H20  pa¬ 
rameterization  for  middle  atmospheric  water  vapor  photo¬ 
chemistry.  CHEM2D-H20  is  based  on  output  from  a  zon- 
ally  averaged  (2-D)  model  with  full  photochemistry  and  ac¬ 
counts  for  the  latitude,  altitude,  and  seasonal  dependences  in 
the  production  and  loss  terms  due  to  CH4  oxidation  and  H2O 
photolysis,  respectively.  The  parameterization  is  valid  up  to 
lower  thermospheric  altitudes  of  ~95-100km. 

A  series  of  10-day  NOGAPS-ALPHA  forecast  model  sim¬ 
ulations  with  parameterized  H2O  photochemistry  were  per¬ 
formed  for  the  June  2007  period  to  evaluate  the  model’s  prog¬ 
nostic  capability  for  middle  atmospheric  H2O.  The  forecast 
model  H2O  mixing  ratios  were  compared  with  NOGAPS- 
ALPHA  analyzed  H2O  fields  based  on  assimilation  of  Aura 
MLS  profile  measurements.  We  find  that  the  CHEM2D-H20 
and  ECMWF  parameterizations  both  perform  comparably  at 
all  altitudes  below  the  0. 1  hPa  level.  Above  this  level,  we  find 
that  CHEM2D-H20  performance  is  dependent  on  the  choice 
of  the  reference  state  mixing  ratio  distribution.  In  the  June 
2007  case  examined  here,  using  CHEM2D-H20  with  a  ref¬ 
erence  state  mixing  ratio  distribution  based  on  the  monthly 
zonal  mean  Aura  MLS  H2O  from  the  June  2005-2008  period 
produced  the  best  overall  agreement  with  the  analyses.  Mid¬ 
dle  atmospheric  water  vapor  mixing  ratios  were  lower  in  the 
early  1990’s  and  exhibited  significant  positive  trends  during 
the  UARS  time  period  (e.g.,  Nedoluha  et  al.,  1998;  Randel 
et  al.,  1999).  Therefore  it  is  not  surprising  that  the  use  of 


the  UARS  H2O  climatology  to  specify  r°  in  EXP1  leads  to 
a  systematic  underprediction  for  the  June  2007  period  ana¬ 
lyzed  here. 

These  preliminary  results  show  that  the  latitude  and  sea¬ 
sonal  dependences  of  the  CHEM2D-H20  parameterization 
can  offer  an  improved  10-day  forecast  of  upper  mesospheric 
water  vapor  compared  to  the  1-D  parameterization  currently 
used  in  some  operational  NWP/DA  systems.  We  are  now 
carrying  out  a  more  formal  evaluation  of  middle  atmospheric 
H2O  prognostic  skill  in  NOGAPS-ALPHA  over  a  longer 
time  period  using  the  fully  coupled  NWP/DA  system,  which 
will  provide  a  more  complete  and  objective  assessment  of 
CHEM2D-H20  performance.  Future  work  will  combine 
NOGAPS-ALPHA  prognostic  H2O  and  temperature  fields 
in  order  to  predict  supersaturation  conditions  that  lead  to  the 
formation  of  polar  mesospheric  clouds.  The  role  of  addi¬ 
tional  physical  processes  such  as  molecular  diffusion,  which 
is  not  currently  accounted  for  in  either  the  CHEM2D  or 
NOGAPS-ALPHA  model,  will  also  be  investigated. 

The  CHEM2D-H20  parameterization  is  freely  available 
for  research  purposes.  For  more  details,  please  contact  the 
lead  author. 
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